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We demonstrate that the dynamics towards and within the Feigenbaum attractor combine to form 
a g-deformed statistical-mechanical construction. The rate at which ensemble trajectories converge 
to the attractor (and to the repellor) is described by a g-entropy obtained from a partition function 
generated by summing distances between neighboring positions of the attractor. The values of the 
g-indices involved are given by the unimodal map universal constants, while the thermodynamic 
structure is closely related to that formerly developed for multifractals. As an essential component 
in our demonstration we expose, in great detail, the features of the dynamics of trajectories that 
either evolve towards the Feigenbaum attractor or are captured by its matching repellor. The 
dynamical properties of the family of periodic superstable cycles in unimodal maps are seen to 
be key ingredients for the comprehension of the discrete scale invariance features present at the 
period-doubling transition to chaos. Elements in our analysis are the following, (i) The preimages 
of the attractor and repellor of each of the supercycles appear entrenched into a fractal hierarchical 
structure of increasing complexity as period doubling develops, (ii) The limiting form of this rank 
structure results in an infinite number of families of well-defined phase-space gaps in the positions 
of the Feigenbaum attractor or of its repellor. (iii) The gaps in each of these families can be 
ordered with decreasing width in accordance with power laws and are seen to appear sequentially 
in the dynamics generated by uniform distributions of initial conditions, (iv) The power law with 
log-periodic modulation associated with the rate of approach of trajectories towards the attractor 
(and to the repellor) is explained in terms of the progression of gap formation, (v) The relationship 
between the law of rate of convergence to the attractor and the inexhaustible hierarchy feature of the 
preimage structure is elucidated, (vi) A "mean field" evaluation of the atypical partition function, a 
thermodynamic interpretation of the time evolution process, and a crossover to ordinary exponential 
statistics are given. We make clear the dynamical origin of the anomalous thermodynamic framework 
existing at the Feigenbaum attractor. 
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I. INTRODUCTION 

A fundamental question in statistical physics is 
whether the structure of ordinary equilibrium statisti- 
cal mechanics falters when its fundamental properties, 
phase-space mixing and ergodicity, breakdown. The 
chaotic dynamics displayed by dissipative nonlinear sys- 
tems, even those of low dimensionality, possesses these 
two crucial conditions, and acts in accordance with a for- 
mal structure analogous to that of canonical statistical 
mechanics, in which thermodynamic concepts meet their 
dynamical counterparts [1]. At the transition between 
chaotic and regular behavior, classically represented by 
the Feigenbaum attractor [1 , the Lyapunov exponent 
vanishes and chaotic dynamics turns critical. Trajecto- 
ries cease to be ergodic and mixing; they retain mem- 
ory of their initial positions and fluctuate according to 
complex deterministic patterns [5]. Under these condi- 
tions it is of interest to check up whether the statistical- 
mechanical structure subsists, and if so, examine if it is 
unchanged or if it has acquired a new form. 

With this purpose in mind, the exploration of possible 



limits of validity of the canonical statistical mechanics, an 
ideal model system is a one-dimensional map at the tran- 
sition between chaotic and regular behavior, represented 
by well-known critical attractors, such as the Feigenbaum 
attractor. So far, recent studies [3j have concentrated on 
the dynamics inside the attractor and have revealed that 
these trajectories obey remarkably rich scaling properties 
not known previously at this level of detail [4,. The re- 
sults are exact and clarify [1] the relationship between 
the original modification [21 |5] of the thermodynamic 
approach to chaotic attractors [51 [71 15] for this type of 
incipiently chaotic attractor, and some aspects of the 
q-deformed statistical mechanical formalism [^ [THl E] • 
The complementary part of the dynamics, that of ad- 
vance on the way to the attractor, has, to our knowl- 
edge, not been analyzed, nor understood, with a simi- 
lar degree of thoroughness. The process of convergence 
of trajectories into the Feigenbaum attractor poses sev- 
eral interesting questions that we attempt to answer here 
and elsewhere |12j based on the comprehensive knowledge 
presented below. Prominent amongst these questions is 
the nature of the connection between the two sets of dy- 
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namical properties, within and outside the attractor. As 
it turns out, these two sets of properties are related to 
each other in a statistical-mechanical manner, i.e. the 
dynamics at the attractor provides the configurations in 
a partition function while the approach to the attractor 
is described by an entropy obtained from it. As we show 
below, this statistical-mechanical property conforms to 
a q-deformation [9J of the ordinary exponential weight 
statistics. 

Trajectories inside the attractor visit positions forming 
oscillating log-periodic patterns of ever increasing am- 
plitude. However, when the trajectories are observed 
only at specified times, positions align according to power 
laws, or g-exponential functions that share the same q- 
index value [H [10]. Further, all such sequences of posi- 
tions can be shifted and seen to collapse into a single one 
by a rescaling operation similar to that observed for cor- 
relations in glassy dynamics, a property known as "ag- 
ing" [ini [13]. The structure found in the dynamics is 
also seen to consist of a family of Mori's q- phase transi- 
tions P^, via which the connection is made between the 
modified thermodynamic approach and the g-statistical 
property of the sensitivity to initial conditions [H [TO] . 
On the other hand, a foretaste of the nature of the dy- 
namics outside the critical attractor can be appreciated 
by considering the dynamics towards the so-called super- 
stable cycles, or supercycles, the family of periodic at- 
tractors with Lyapunov exponents that diverge towards 
minus infinity. This infinite family of attractors has as 
accumulation point the transition to chaos, which for the 
period-doubling route is the Feigenbaum attractor. As 
described below, the basins of attraction for the differ- 
ent positions of the cycles develop fractal boundaries of 
increasing complexity as the period-doubling structure 
advances towards the transition to chaos. The fractal 
boundaries, formed by the preimages of the repellor, dis- 
play hierarchical structures organized according to ex- 
ponential clusterings that manifest in the dynamics as 
sensitivity to the final state and transient chaos. The 
hierarchical arrangement expands as the period of the 
supercycle increases. 

To observe the general procedure followed by trajecto- 
ries to reach the attractors, and their complementary re- 
pellors, we consider an ensemble of uniformly distributed 
initial conditions Xq spanning the entire phase space in- 
terval. We find that this is a highly structured process 
encoded in sequences of positions shared by as many tra- 
jectories with different xq. Clearly, there is always a nat- 
ural dynamical ordering in the xq as any trajectory of 
length t contains consecutive positions of other trajecto- 
ries of lengths t ~ I, t — 2, etc. with initial conditions 
Xq, Xq, etc. that are images under repeated map itera- 
tions of Xq- In the case of the Feigenbaum attractor the 
initial conditions form two sets, dense in each other, of 
preimages of each the attractor and the repellor. There 
is an infinite-level structure within these sets that, as 
we shall see, is reflected by the infinite number of fami- 
lies of phase-space gaps that complement the multifrac- 



tal layout of both attractor and repellor. These families 
of gaps appear sequentially in the dynamics, beginning 
with the largest and followed by other sets consisting of 
continually increasing elements with decreasing widths. 
The number of gaps in each set of comparable widths 
increases as 2*^, fc = 0, 1, . . . and their widths can be or- 
dered according to power laws of the form a~'', where 
a is Feigenbaum's universal constant a ~ 2.5091. We 
call k the order of the gap set. Furthermore, by consid- 
ering a fine partition of phase-space, we determine the 
overall rate of approach of trajectories towards the at- 
tractor (and to the repellor). This rate is measured by 
the fraction of bins W{t) still occupied by trajectories at 
time t |14j . The power law with log-periodic modulation 
displayed by W{t) [H] is explained in terms of the pro- 
gression of gap formation, and its self-similar features are 
seen to originate in the unlimited hierarchy feature of the 
preimage structure. 



Before proceeding to give details in the following sec- 
tions of the aforementioned dynamics, we recall basic 
features of the bifurcation forks that form the period- 
doubling cascade sequence in unimodal maps, epitomized 
by the logistic map /^(x) = 1 — fix^, — 1 < x < 1, 
< ^ < 2 [I] |T5|. The superstable periodic orbits of 
lengths 2^, N = 1,2,3,..., are located along the bifurca- 
tion forks, i.e. the control parameter value /i = /J^ < /ioo 
for the superstable 2^-attractor is that for which the or- 
bit of period 2^ contains the point x = 0, where fioo — 
1.401155189. . . is the value of at the period-doubling 
accumulation point. The positions (or phases) of the 2^- 
attractor are given by Xm = /p™^(0), m = 1, 2, ... , 2^. 
Notice that infinitely many other sequences of super- 
stable attractors appear at the period-doubling cascades 
within the windows of periodic attractors for values of 
A* > Moo- Associated with the 2^-attractor at /i = JIj^ 
there is a (2^ — l)-repellor consisting of 2^ — 1 posi- 



tions Urn, m = 1,2, 



unstable solutions 
fc = 0,l,2, ...,7V- 



1. These positions are the 



1. The first, fc = 0, originates at the 
initial period-doubling bifurcation, the next two, fc = 1, 
start at the second bifurcation, and so on, with the last 
group of 2^~^, fc = — 1, stemming from the — 1 
bifurcation. We find it useful to order the repellor posi- 
tions, or simply, repellors, present a.t fj, — JLj^, according 
to a hierarchy or tree, the "oldest" with fc = up to 
the most "recent" ones with fc = A^ — 1. The repel- 
lors' order is given by the value of fc. Finally, we de- 
fine the preimage x'^^ of order j of position x to satisfy 
X = h^^\x'^^'>) where h'^^\x) is the j-ih composition of 
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the map h{x) = /A^ We have omitted reference 

to the 2^-cycle in x^^^ to simplify the notation. The 
interval lengths or diameters diq^m are measured when 
considering the superstable periodic orbits of lengths 2^. 
The di^^m are defined (here) as the (positive) distances of 



the elements Xm, m = 0, 1, 2, 



1, to their nearest 
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neighbors i.e. 



(1) 



For large N, dN,o/dN+i.o — ct- We present explicit re- 
sults for the logistic map, which has a quadratic max- 
imum, but the results are easily extended to unimodal 
maps with general nonlinearity z > 1. 

Central to our discussion is the following broad prop- 
erty: Time evolution at fi^o from t = up to i — > cx3 
traces the period-doubling cascade progression from fi = 
up to fioo- Not only is there a close resemblance be- 
tween the two developments but also quantitative agree- 
ment. For instance, the trajectory inside the Feigenbaum 
attractor with initial condition xq = 0, the 2°°-supercycle 
orbit, takes positions Xt such that the distances between 
appropriate pairs of them reproduce the diameters dN^m 
defined from the supercycle orbits with /I^ < fioo- See 
Fig. [l] where the absolute value of positions and loga- 
rithmic scales are used to illustrate the equivalence. This 
property has been key to obtain rigorous results for the 
sensitivity to initial conditions for the Feigenbaum at- 
tractor iSillDj. 
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FIG. 1: Left panel: Absolute value of attractor positions 
for the logistic map /^(x) in logarithmic scale as a function 
of the logarithm of the control parameter difference /ioo— A*- 
Right panel: Absolute value of trajectory positions for the 
logistic map ffxix) at fi^o with initial condition 2:0= in log- 
arithmic scale as a function of the logarithm of time t, also 
shown by the numbers close to the circles. The arrows in- 
dicate the equivalence between the diameters 5m in the left 
panel, and position differences Sn with respect to 2;o= in 
the right panel. 

The layout of the rest of this paper is the following. In 
Sec. [njwe present the dynamical properties of the fam- 
ily of supercycles, describing their preimage structure, 
final state sensitivity and transient chaos. Details of the 
superstrong insensitivity to initial conditions displayed 



by these attractors are given as an Appendix. In Sec. 
|III| we make use of the results of the previous section 
to describe the dynamical properties of approach to the 
Feigenbaum attractor. We provide details of its preim- 
age structure, the sequential opening of phase-space gaps, 
and the scaling for the trajectories' rate of convergence to 
the attractor (and repellor). In Sec. IV we explain the 
aforesaid statistical-mechanical structure lying beneath 
the dynamics of an ensemble of trajectories en route to 
the Feigenbaum attractor (and repellor). In Sec. [v| we 
summarize our results. 



II. HIERARCHICAL PROPERTIES IN THE 
DYNAMICS OF SUPERCYCLE ATTRACTORS 

To obtain dynamical properties with previously un- 
stated detail we determined the organization of the en- 
tire set of trajectories as generated by all possible initial 
conditions. We find that the paths taken by the full set 
of trajectories in their way to the supercycle attractors 
(or to their complementary repellors) are far from un- 
structured. The preimages of the attractor of period 2''^, 
TV = 1, 2, 3, . . . are distributed into different basins of at- 
traction, one for each of the 2^ phases (positions) that 
compose the cycle. When N > 2 these basins are sep- 
arated by fractal boundaries whose complexity increases 
with increasing N . The boundaries consist of the preim- 
ages of the corresponding repellor and their positions 
cluster around the 2^ — 1 repellor positions according to 
an exponential law. As N increases the structure of the 
basin boundaries becomes more involved. That is, the 
boundaries for the 2''^ cycle develop new features around 
those of the previous 2^^^ cycle boundaries, with the 
outcome that a hierarchical structure arises, leading to 
embedded clusters of clusters of boundary positions, and 
so forth. 

The dynamics associated with families of trajectories 
always displays a distinctively concerted order that re- 
flects the repellor preimage boundary structure of the 
basins of attraction. That is, each trajectory has an ini- 
tial condition that is identified as an attractor (or repel- 
lor) preimage of a given order, and this trajectory neces- 
sarily follows the steps of other trajectories with initial 
conditions of lower preimage order belonging to a given 
chain or pathway to the attractor (or repellor). This 
feature gives rise to transient chaotic behavior different 
from that observed at the last stage of approach to the 
attractor. When the period 2^ of the cycle increases 
the dynamics becomes more involved with increasingly 
more complex stages that reflect the preimage hierarchi- 
cal structure. As a flnal point, shown in the appendix, 
in the closing part of the last leg of the trajectories an 
ultra-rapid convergence to the attractor is observed, with 
a sensitivity to initial conditions that decreases as an ex- 
ponential of an exponential in time. In relation to this 
we flnd that there is a functional composition renormal- 
ization group (RG) flxed-point map associated with the 
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supercycle attractor, and this can be expressed in closed 
form by the same kind of q -exponential function found 
for both the pitchfork and tangent bifurcation attractors 
[ini HZ] ) like that originally derived by Hu and Rudnick 
for the latter case [18] , 



A. Preimage structure of supercycle attractors 

The core source of our description of the dynamics to- 
wards the supercycle attractors is a measure of the rel- 
ative "time of flight" tf{xo) for a trajectory with initial 
condition xq to reach the attractor. The function tf^xo) 
is obtained for an ensemble representative of all initial 
conditions — 1 < xq < 1. This comprehensive informa- 
tion is determined through the numerical realization of 
every trajectory, up to a small cutoff e > at its final 
stage. The cutoff e considers a position Xf < Xm ± e to 
be effectively the attractor phase Xm- This, of course, 
introduces an approximation to the real time of flight, 
which can be arbitrarily large for those xq close to a re- 
pellor position or close to any of its infinitely many 
preimages, xin , j — 1,2,..., N > 2. In such cases 
the finite time tf(xo;s) can be seen to diverge tf ^ oo 
as Xq ^ Um and e ^ 0. As a simple illustration, in 
Fig. |2]we show the time of fiight tf{xo) for the period 
2 supercycle at ^ = with e = 10~^ together with 

(2) 

the twice-composed map fi.^'[x) and a few representa- 
tive trajectories. We observe two peaks in tf{xo;e) at 
2/1 = (-1 + VI +4a*i)/2/Ii ~ 0.6180340..., the fixed- 
point repellor and at its (only) preimage x^i ^ — —yi, 
where j/i = f'^^{x^i'^). Clearly, there are two basins of 
attraction each for the two positions or phases, xi — Q 
and X2 = 1, of the attractor. For the former it is the 
interval x^^'' < Xq < yi, whereas for the second one it 

consists of the intervals — 1 < xo < and yi < x^ <\. 
The multiple-step structure of tf{xo), of four time units 
each step, reflects the occurrence of intervals of initial 
conditions with common attractor phase preimage order 
k. 

For the next supercycle — period 4 — the preimage 
structure turns out to be a good deal more involved than 
the straightforward structure for /Zj^. In Fig. [3] we show 
the times of flight tf{xo) for the N = 2 supercycle at 

/i = ^2 ^itli £ = 10^^, the map f^^{x) is superposed as 
a reference to indicate the four phases of the attractor (at 
= 0, 2)2 = 1, xa ~ -0.3107 . . ., and xa ~ 0.8734 . . .) 
and the three repellor positions (at yi ~ 0.5716635..., 
?/2 ^ 0.952771..., and 2/3 ~ -0.189822...) In Fig. [s] 
there are also shown four trajectories each of which ter- 
minates at a different attractor phase. We observe a pro- 
liferation of peaks and valleys in tf^xo), actually, an infi- 
nite number of them, that cluster around the repellor at 
2/1 ~ 0.5716635 . . . and also at its preimage at x^^^ — —yi 
(these are the positions at = /Jj of the "old" repellor 
and its preimage in the previous = 1 case). Notice 




FIG. 2: Left axis: The twice iterated map {x), with /ij^= 1 
(gray line). Right axis: Time of flight tj{x), the number of 
iterations necessary for a trajectory with initial condition at 
X to reach an attractor position. The values of x near the 
high spikes correspond to initial conditions very close to the 
repellor and its preimage. We present three example trajec- 
tories (and the y ~ x line as an aid to visualize them): The 
dotted line shows a trajectory that starts at the attractor po- 
sition a; = and remains there. The solid line is a trajectory 
starting near the repellor at yi, and after a large number of 
iterations reaches the attractor position x = 1. Finally, the 
dash-dotted line is an orbit starting at a; = 0.5 that in just a 
few iterations reaches a; = 0. 



that the steps in the valleys of tf{xo) are now eight time 
units each. The nature of the clustering of peaks (re- 
pellor phase preimages) and the bases of the valleys (at- 
tractor phase preimages) is revealed in Fig. |4] where we 
plot on a logarithmic scale for the variables ±(a; — 2/1). 
There is an exponential clustering of the preimage struc- 
ture around both the old repellor and its preimage. This 
scaling property is corroborated in Fig. |5]from which we 
obtain x-2/1 ~ 7.5x 10"^ exp(0.80 0, where / = 1,2, 3, . . . 
is a label for consecutive repellor preimages. 

A comparable leap in the complexity of the preim- 
age structure is observed for the following — period 8 — 
supercycle. In Fig. |6] we show tf{xo) for the N — 3 
supercycle at /i = pa with e — 10~^, together with the 
map fj?^(^x) placed as reference to facilitate the identi- 
fication of the locations of the eight phases of the at- 
tractor, xi to xg, and the seven repellor positions, 2/1 
to 2/7. In addition to a huge proliferation of peaks and 
valleys in tf(xo), we observe now the development of 
clusters of clusters of peaks centered around the repel- 
lor at 2/1 — 0.5626447. . . and its preimage at xj^^ — —yi 
(these are now the positions of the original repellor and 
its preimage for the N = 1 case when fj. — JL^). The 
steps in the valleys of i/(xo) have become 16 time units 
each. Similarly to the clustering of peaks (repellor phase 
preimages) and valleys (attractor phase preimages) for 
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FIG. 3: Same as Fig. M but for f^Jix), JI^ ~ 1.31070264. 
Two orbits (solid bold and dashed bold lines) start very near 
each other and by the position yi of the old repellor, although 
indistinguishable at the beginning they take very different 
subsequent paths to reach the attractor positions at a; = 1 
and X ~ 0.8734. See bottom right panel. Another two or- 
bits (in dotted and solid lines) start at an attractor position 
X = 1 and at a repellor preimage position, respectively. Fi- 
nally, one more orbit (dashed line) starts at an intermediate 
initial condition and reaches very quickly the attractor posi- 
tion at a; ~ —0.310703. See bottom left panel. 



the previous supercycle at /Ij, the spacing arrangement 
of the new clusters of clusters of peaks is determined in 
Fig. It] where we plot tf in a. logarithmic scale for the 
variables zt{x — yi) . In parallel to the previous cycle an 
exponential clustering of clusters of the preimage struc- 
ture is found around both the old repellor and its preim- 
age. This scaling property is quantified in Fig. [8] from 
which we obtain a; — j/i ~ 8.8 x 10"^ exp(0.84/), where 
I — 1,2,3, . . . counts consecutive clusters. 

An investigation of the preimage structure for the next 
= 4 supercycle dX ji — ~p.^ leads to another substan- 
tial increment in the complications of the structure of 
the preimages but with such density that is cumbersome 
to describe here. Nevertheless it is clear that the main 
characteristic in the dynamics is the development of a hi- 
erarchical organization of the preimage structure as the 
period 2^ of the supercycles increases. 



B. Final state sensitivity and transient chaos 

With the knowledge gained about the features dis- 
played by the times of flight t f (xq) for the first few super- 
cycles it is possible to determine and understand how the 
leading properties in the dynamics of approach to these 
attractors arise. The information contained in tf{xQ) can 
be used to demonstrate in detail how the concepts of final 




>'j-X X-y^ 

FIG. 4: Top panel: Time of flight tf{x) for AT = 2 (as in Fig. 
3), the black lines correspond to initial conditions that termi- 
nate at the attractor positions a; = and x ~ —0.310703, and 
the gray lines to trajectories ending at a; = 1 and x ~ 0.8734. 
Right (left) bottom panel: Same as top panel, but plotted 
against the logarithm of x — [y\~x). It is evident that the 
peaks are arranged exponentially around the old repellor po- 
sition j/i, i.e., they appear equidistant on a logarithmic scale. 
See Fig. [5] 
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FIG. 5: Corroboration of preimage exponential clustering 
around the repellor position j/i ~ 0.571663 when N = 2 and 
'p2 — 1.31070264. The variable I labels consecutive equidistant 
peaks in the inset. The peaks correspond to preimages of the 
repellor. See text. 



state sensitivity [12] — due to attractor multiplicity — and 
transient chaos [1] — prevalent in the presence of repel- 
lors that coexist with periodic attractors — are realized 
in a given dynamics. Final state sensitivity is the con- 
sequence of fractal boundaries separating coexisting at- 
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FIG. 6: Same as Fig. 3 but for f^l'{x), Jl^~ 1.38154748. Here 
2/1= 0.56264475. 




FIG. 7: Same as Fig. |4] but for iV = 3. The black lines 
correspond to initial conditions that terminate at any of the 
four attractor positions close or equal to a; = 0, and the gray 
lines to trajectories ending at any of the other four attractor 
positions close or equal to x = 1. As the bottom panels show, 
on a logarithmic scale, in this case there are (infinitely) many 
clusters of peaks (repellor preimages) equidistant from each 
other. 



tractors. In our case there is always a single attractor but 
its positions or phases play an equivalent role [20 . Tran- 
sient chaos [1 is due to fast separation in time of nearby 
trajectories by the action of a repellor and results in a 
sensitivity to initial conditions that grows exponentially 
up to a crossover time after which decay sets in. We 
describe how both properties result from an extremely 
ordered flow of trajectories towards the attractor. This 
order is imprinted by the preimage structure described 
above. 




50 100 150 200 250 300 
t 

FIG. 9: Top panel: Detail of the map fj^J{x) and a tra- 
jectory in the proximity of the attractor position located at 
X ~ 0.850780. The trajectory originated very close to the re- 
pellor position a,t X — 0.819378. Bottom panel: The same 
trajectory of the eight-iterated map as a function of time t. 
Inset: Corroboration of the exponential nature of the trajec- 
tory after it leaves the repellor and before the final approach 
to the attractor position. 

For the simplest supercycle at /i = T^i there is trivial 
final state sensitivity as the boundary between the two 
basins of the phases, xi = and X2 = 1, consists only of 
the two positions yi ~ 0.5716635 . . . and xj^'' = —yi- See 
Fig. [2j Consider the length 5 of a small interval around a 
given value of containing either yi or x\^\ when 5 — *■ 
any uncertainty as to the final phase of the trajectory 
disappears. It is also simple to verify that when xq is 
close to 2/1 or x^^'' the resulting trajectories increase their 
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separation at initial and intermediate times displaying 
transient chaos in a straightforward fashion. In Fig. [9] we 
show the same type of transitory exponential sensitivity 
to initial conditions for a trajectory ai fi = JI^ after it 
reaches a repellor position in the final journey towards 
the period 8 attractor. This behavior is common to all 
periodic attractors when trajectories come near a repellor 
at the final leg of their journey. 
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FIG. 10: Time of flight tf{x) for iV = 2 (equivalent to Fig. 
|4|. Different types of lines are used for different sub-basins of 
attraction that form the fractal boundary between the attrac- 
tor positions. The black solid line corresponds to the attrac- 
tor position X = 0, the gray solid line to x = 1, the dashed 
gray line to a: ~ —0.310703, and the dashed black line to 
X ~ 0.873470. The sub-basins are separated in an alternating 
fashion by repellor preimages (peaks) and there are just two 
types of sub-basins depending on |a::| being larger or smaller 
than yi. 

For fi ~ J12 there are more remarkable properties aris- 
ing from the more complex preimage structure. There 
is a concerted migration of initial conditions seeping 
through the boundaries between the four basins of attrac- 
tion of the phases. These boundaries, shown in Fig. [TO] 
form a fractal network of interlaced sub-basins separated 
from each other by two preimages of different repellor 
phases and have at their bottom a preimage of an attrac- 
tor phase. Trajectories on one of these sub-basins move 
to the nearest sub-basin of its type (next-nearest neigh- 
bor in actual distance in Fig 



map fj^^ (four time steps for the original map /jr^ ) 
movement is always away from the center of the cluster 
at the old repellor position j/i or at its preimage x']^^ (lo- 
cated at the steepest slope infiection points of fj^^ shown 
in Fig. [3]). Once a trajectory is out of the cluster (con- 
tained between the maxima and minima of fi*-* next to 
the mentioned inflection points) it proceeds to the basin 
of attraction of an attractor phase (separated from the 
cluster by the inflection points with gentler slope of fj^^^ 



10 ) at each iteration of the 
The 



in Fig. |3| where its final stage takes place. When we 
consider a large ensemble of initial positions, distributed 
uniformly along all phase-space, the common journey to- 
wards the attractor displays an exceedingly ordered pat- 
tern. Each initial position xq within either of the two 
clusters of sub-basins is a preimage of a given order k of 
a position in the main basin of attraction. Each iteration 
of fj^^ reduces the order of the preimage from fc to fc — 1 , 

and the new position Xq = /l'^''(xo) replaces the initial 
position Xq (a preimage of order k — I) of another trajec- 
tory that under the same time step has migrated to the 
initial position Xq (a preimage of order fc — 2) of another 
trajectory, and so on. 




FIG. 11: Evolution of an orbit starting very close to a 
peak (repellor preimage) inside a cluster when N = 2, JI2 — 
1.31070264 and yi= 0.571663. (See Fig. [I]) The trajectory 
describes an exponential in time moving away from yi. In 
the figure, positions of the trajectory correspond to the black 
solid dots, and iterations correspond to the number associated 
with each dot. 

It is clear that the dynamics at ^ = ^2 displays sen- 
sitivity to the final state when the initial condition xq is 
located near the core of any of the two clusters of sub- 
basins at yi and x^^^ that form the boundary of the at- 
tractor phases. Any uncertainty on the location of Xq 
when arbitrarily close to these positions implies uncer- 
tainty about the final phase of a trajectory. There is also 
transient chaotic behavior associated to the migration of 
trajectories out of the cluster as a result of the organized 
preimage resettlement mentioned above. Indeed, the ex- 
ponential disposition of repellor preimages shown in Fig. 
[5] is actually a realization of two trajectories with ini- 
tial conditions in consecutive peaks of the cluster struc- 
ture. Therefore, the exponential expression given in the 
previous section in relation to Fig. [5] can be rewritten 
as the expression for a trajectory xg ~ exp(Ae//6'), 
with Xg — X ~- ui, Xq — 7.5 X 10~^, and Ag// = 6.4, 
where 0=1,2,3,... Straightforward differentiation of xg 
with respect to xq yields an exponential sensitivity to 
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initial conditions with positive effective Lyapunov coeffi- 
cient Ae//. See Fig. 



11 



As can be anticipated, the dynamics of approach to the 
next supercycle a,t ^ = JL^ can be explained by enlarging 
the description presented above for fi = ]l2 with the addi- 
tional features of its preimage structure already detailed 
in the preceding section. As in the previous case, trajec- 
tories with initial conditions xq located inside a cluster of 
sub-basins of the attractor phases will proceed to move 
out of it in the systematic manner described for the only 
two isolated clusters present when jj. = 'p2. However, now 
there is an infinite number of such clusters arranged into 
two bunches that group exponentially around the old re- 
pellor position yi and around its preimage x'^^K See Figs. 
[6}|8] Once such trajectories leave the cluster under con- 
sideration they enter into a neighboring cluster, and so 
forth, so that the trajectories advance out of these fractal 
boundaries through the prolonged process of migration 
out of the cluster of clusters before they proceed to the 
basins of the attraction of the eight phases of this cycle. 
In Fig. [12] we show one such trajectory in consecutive 
times i = 1, 2, 3, . . . for the original map and also in mul- 
tiples of time i = 2^, 22^, 32^, .... The logarithmic scale 
of the figure makes evident the retardation of each stage 
in the process. As when ^ = 'p2, it is clear that in the 
approach to the fx = Ji^ attractor there is sensitivity to 
the final state and transitory chaotic sensitivity to initial 
conditions. Again, the exponential expression, given in 
the previous section associated with the preimage struc- 
ture of clusters of clusters of sub-basins, shown in Fig. 
[S] can be interpreted as the expression of a trajectory of 
the form xg ~ xq exp(Ae//0). Differentiation of xe with 
respect to xq yields again an exponential sensitivity to 
initial conditions with positive effective Lyapunov coeffi- 
cient Ae//. 

When the period 2^ of the cycles increases we observe 
that the main characteristic in the dynamics is the de- 
velopment of a hierarchical organization in the flow of an 
ensemble of trajectories out of an increasingly more com- 
plex disposition of the preimages of the attractor phases. 



III. DYNAMICAL PROPERTIES OF 
APPROACH TO THE FEIGENBAUM 
ATTRACTOR 

A convenient way to visualize how the preimages for 
the Feigenbaum attractor and repellor are distributed 
and organized is to consider the simpler arrangements 
for the preimages of the supercycles' attractors and re- 
pellors of small periods 2^, N = 1,2,... As we have 
seen, when the period 2^ increases the preimage struc- 
tures for the attractor and repellor become more and 
more involved, with the appearance of new features made 
up of an infinite repetition of building blocks. Each of 
these new blocks is equivalent to the more dense struc- 
tures present in the previous 2^~^ case. In addition all 
other structures in the earlier 2^~^, 2^ cases are still 
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FIG. 12: A trajectory for Ji^~ 1.38154748. Top panel: The 
circles are positions for consecutive times in the iterations of 
the map /jj^ (x). The orbit starts very close to the old (period 
1) repellor yi, then moves close to a period 2 and subsequently 
to a period 4 repellor, until finally it arrives at the period 8 
attractor. Bottom panel: Selection of a time subsequence 
of multiples of 2^ shows the same trajectory as an evolution 
towards one particular attractor position or phase. 



present. Thus a progressively more elaborate organiza- 
tion of preimages is built upon as N increases, so that the 
preimage layout for the Feigenbaum attractor and repel- 
lor is obtained as the limiting form of the rank structure 
of the fractal boundaries between the finite period at- 
tractor position basins. The fractal boundaries consist of 
sub-basins of preimages for the attractor positions sep- 
arated by preimages of the repellor positions. As N in- 
creases the sizes of these sub-basins decrease while their 
numbers increase and the fractal boundaries cover a pro- 
gressively larger part of total phase-space. See Figs. [T3| 
andd 

Interestingly, the sizes of all boundary sub-basins van- 
ish in the limit N —f 00, and the preimages of both 
attractor and repellor positions become two sets — with 
dimension equal to the dimension of phase space — dense 
in each other. In the limit N 00 there is an attractor 
preimage between any two repellor preimages and the 
other way round. (The attractor and repellor are two 
multifractal sets with dimension d/ — 0.538... [E].) To 
visualize this limiting situation consider that the posi- 
tions for the repellors and their first preimages of the 
2^-th supercycle appear located at the infiection points 

of fjj-^ (x) , and it is in the close vicinity of them that the 
mentioned fractal boundaries form. To illustrate how the 
sets of preimage structures for the Feigenbaum attractor 
and repellor develop we plot in Fig. [13] the absolute value 

of In 



41 Vrf- 



for N — 1,2, ...,4 vs x. The maxima in 



(2 ) 

this curve correspond to the inflection points of /^^ (x) 
at which the repellor positions or their first preimages 
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are located. As shown in Fig. 13 when N increases the 
number of maxima proliferate at a rate faster than 2^. 




FIG. 13: Absolute value of df'^^\x) /dx, for = 1, 2, 3 and 
4, on a logarithmic scale as a function of x in the interval 
< a; < 1. The proliferation of maxima conveys the develop- 
ment of the hierarchical structure of repellor preimages. See 
text. 



A. Sequential opening of phase-space gaps 

One way that the preimage structure described above 
is manifest in the dynamics is via the successive formation 
of phase-space gaps that ultimately give rise to the at- 
tractor and repellor multifractal sets. In order to observe 
explicitly this process we consider an ensamble of initial 
conditions spread out uniformly across the interval 
— 1 < < 1 and keep track of their positions at sub- 
sequent times. In Figs. 



14 



16 we illustrate the outcome 
for the supercycles of periods 2^, 2'^, and 2^, respectively, 
where we have plotted the time evolution of an ensemble 
composed of 10'' trajectories. In the left panel of each 
figure we show the absolute value of the positions \xt \ vs 
time t, while, for comparison purposes, in the right panel 

we show the absolute value of both vs /^^ {x) and 



to facilitate identification of the attractor 

and repellor positions. The labels fc = 0, 1, 2, . . . indicate 
the order of the gap set (or equivalently the order of the 
repellor generation set). In Fig. [14] (with /i = o'^g 
observes a large gap opening first that contains the old 
repellor (fc = 0) in its middle region and two smaller gaps 
opening afterward that contain the two repellors of sec- 
ond generation {k — 1) once more around the middle of 
them. In Fig. 15 (with — /Ig) we initially observe the 
opening of a primary and the two secondary gaps as in the 
previous /i = /i2 case, but subsequently four new smaller 
gaps open each around the third generation of repellor 
positions (fc = 2). In Fig. 16 (with ^ = ^4) we observe 



the same development as before; however, at longer times 
eight additional and yet smaller gaps emerge around each 
fourth generation of repellor positions (fc = 3). Naturally, 
this process continues indefinitely as A'^ — *■ cxd and illus- 
trates the property mentioned before for /Xqo, that the 
time evolution at fixed control parameter values resem- 
bles progression from /i = up to, in this section, /I^. 
It is evident in all Figs. [T4}|T6| that the closer the ini- 
tial conditions Xq are to the repellor positions the longer 
times it takes for the resultant trajectories to clear the 
gap regions. This intuitively evident feature is essentially 
linked to the knowledge we have gained about the fractal 
boundaries of the preimage structure, and the observable 
"bent over" portions of these distinct trajectories in the 
figures correspond to their passage across the boundaries. 
(Since the ensemble used in the numerical experiments is 
finite there appear only a few such trajectories in Figs. 
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FIG. 14: Phase-space gap formation for /i =Jl2- Left panel: 
time evolution of a uniform ensemble of 10* trajectories as 
a function of |x| (black areas and open circles). The values 
of the index fc label the order of the gap set. Right panel: 

Rotated plots of /p*^ (a;) (gray) and |d/p^' (a;)/cixj (black) vs 

jx| as guides for the identification of attractor and repellor 
positions. 



To facilitate a visual comparison between the process 
of gap formation at [i^a and the dynamics inside the 
Feigenbaum attractor — as illustrated by the trajectory 
in Fig. [1] (right panel) — we plot in Fig. 17 the time evo- 
lution of the same ensemble composed of 10"' trajectories 
with II = fiao- This time we use logarithmic scales for 
both \xt\ and t and then superpose on the evolution of 
the ensemble the positions for the trajectory starting at 
xq = 0. It is clear from this figure that the larger gaps 
that form consecutively all have the same width in the 
logarithmic scale of the plot and therefore their actual 
widths decrease as a power law, the same power law fol- 
lowed, for instance, by the position sequence Xt — , 
t — 2^, A^ = 0, 1, 2, . . . for the trajectory inside the at- 
tractor starting at a;o = 0. This set of gaps develop in 
time beginning with the largest one containing the fc = 
repellor, then followed by a second gap, one of a set of 
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two gaps associated with the fc = 1 repellor, next a third 
gap, one gap of a set of four gaps associated with the 
fc = 2 repellor, and so forth. The locations of this spe- 
cific family of consecutive gaps advance monotonically 
towards the sparsest region of the multifractal attrac- 
tor located at a: = 0. The remaining gaps formed at 
each stage converge, of course, to locations near other 
regions of the multifractal, but are not easily seen in Fig. 
17 because of the specific way in which this has been 
plotted (and because of the scale used). 
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In Fig. 

plot the same data differently, with the variable In ]xf re- 
placed by In 1 1 — a; I , where now another specific family of 
gaps, one for each value of fc = 0,1,2,.. ., appears, all 
with the same width on the logarithmic scale; their ac- 
tual widths decrease now as q;~^^, N = 0,1,2,.... The 
locations of this second family of consecutive gaps ad- 
vance monotonically towards the most crowded region of 
the multifractal attractor located a,t x = 1. The time 
necessary for the formation of successive gaps of order 
fc — 0,1,2,... increases as 2'^ because the duration of 
equivalent movements of the trajectories across the corre- 
sponding preimage structures involve the 2'^-th composed 

function Ux). 



k={ k=2 



0.2 








1 k=Q 








/ 


ml k=^ 


k=2 



0.8- 
0.6 



-0.2- 



8!&i-'i'i'i''i'-il 



100 200 
t 



300 



\ ' ' ' I 

_\ ^ : 


' 1 'r- 
> W- 

Jl 

/ - 




} 

\ i\ 
Jf 1 

^1 


- , :rr^ 


1 

,11 



-0.5 



0.5 




100 200 
t 



300 °-l 



-0.5 



0.5 



FIG. 16: Phase-space gap formation for ^ — JL^. Left panel: 
time evolution of a uniform ensemble of 10* trajectories as 
a function of |a;| (black areas and open circles). The values 
of the index k label the order of the gap set. Right panel: 



Rotated plots of fj^^^x) (gray) and jd/A^"' (x) /da; j (black) vs 

|x| as guides for the identification of attractor and repellor 
positions. 
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FIG. 15: 



Phase-space gap formation for ^ =fJ'3- Left panel: 

-,4 



a function of \x\ (black areas and open circles). The values 
of the index k label the order of the gap set. Right panel: 



Rotated plots of /i^ (a;) (gray) and jd/l^ (x) /da; | (black) vs 

|a;| as guides for the identification of attractor and repellor 
positions. 



FIG. 17: Phase-space gap formation for /x = ^ao- The black 
dots correspond to time evolution of a uniform ensemble of 
10* trajectories as a function of \x\ vs t, both on logarithmic 
scales. The open circles are the positions, labeled by the 
times at which they are reached, for the trajectory inside the 
Feigenbaum attractor with initial condition a;o= 0, the same 
as the right panel in Fig. [l] 



B. Scaling for the rates of convergence to the 
attractor and repellor 

There is [l4] an all-inclusive and uncomplicated way 
to measure the rate of convergence of an ensemble of 
trajectories to the attractor (and to the repellor) that 
consists of a single time-dependent quantity. A parti- 
tion of phase-space is made of Ni, equally sized boxes or 
bins and a uniform distribution, of initial conditions 
placed along the interval — 1 < a; < 1, is considered again. 



The number r of trajectories per box is r = Nc/N^. The 
quantity of interest is the number of boxes W(t) that 
contain trajectories at time t. This is shown in Fig. [19] 
on logarithmic scales for the first five supercycles of peri- 
ods 2^ to 2^ where we can observe the following features: 
In all cases W{t) shows a similar initial nearly constant 
plateau [W{t) ~ A, 1 < h < to, h = 0(1)] and a fi- 
nal well-defined decay to zero. As it can be observed in 
the left panel of Fig. [l9] the duration of the final decay 
grows (approximately) proportionally to the period 2^ 
of the supercycle. There is an intermediate slow decay 
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FIG. 18: Same as Fig 

|1 — a;|. Notice the change in slope with respect to Fig, 
the opening of gaps and in the layout of the positions for the 
trajectory inside the attractor. 



of W{t) that develops as N increases with duration also 
(just about) proportional to 2^. For the shortest period 
2^, there is no intermediate feature in W{t); this appears 
first for period 2^ as a single dip and expands with one 
undulation every time N increases by one unit. The ex- 
panding intermediate regime exhibits the development of 
a power law decay with the logarithmic oscillations char- 
acteristic of discrete scale invariance [H]. Clearly, the 
manifestation of discrete invariance is expected to be as- 
sociated with the period-doubling cascade. In the right 
panel of Fig. [T9]wc show a superposition of the five curves 
in Fig. [T9](lcft panel) obtained via rescaling of both W{t) 
and t for each curve according to repeated scale factors. 

The limiting form W(t) for iV — > cx) is shown in the 
left panel of Fig. |20] for various values of r while in its 
right panel we show, for r = 100, a scale amplification of 
W{t) with the same factors employed in Fig. [l9 for the 
supercycles with small periods. The behavior of W{t) at 
/ioo was originally presented in Ref . [14j , where the power 
law exponent ip and the logarithmic oscillation parameter 
A in 



W{t) ~ Ah 



In T 
hTX 



t-tn 



(2) 



were obtained numerically with a precision that corre- 
sponds to r = 10. In Eq. ^ h{x) is a periodic function 
with h{\) — 1 and A is the scaling factor between the 
periods of two consecutive oscillations. More recently, in 
Ref. [22 , it was pointed out that numerical estimates of 
W{t) are subject to large finite-size corrections, and, also, 
that W{t) should scale with the intervals in the triadic 
cantor set construction of the Feigenbaum attractor (25], 
from which the value for = 0.800138194 is reported. 
The values for the rescaling factors in our Figs. [19] and 
|20]suffer from these large finite-size effects due to the rel- 
atively small values of r used in the calculations. This is 



evident since the time scaling factor obtained from these 
data differs by 10% from the exact value of A = 2 implied 
by the discrete scale invariance property created by the 
period-doubling cascade. In Fig. [21] we show the rate 
W{t) and the superposition of repeated amplifications of 
itself (as in the right panel of Fig. 20 1 for increasing val- 
ues of Nc- We find that the scaling factor A converges to 
its limit A = 2. 




FIG. 19: Left panel: Rate W{t), divided by the number 
of boxes Nb employed, of approach to the attractor for the 
supercycles of periods 2^, JV = 1, 2, 3, 4 and 5 on logarithmic 
scales. The expression shown corresponds to the power-law 
decay of the developing logarithmic oscillations. Right panel: 
Superposition of the five curves for W{t) in the left panel via 
n-times repeated rescaling factors shown for the horizontal x 
and vertical y axes. 



We are now in a position to appreciate the dynamical 
mechanism at work behind the features of the decay rate 
W{t). From our previous discussion we know that, every 
time the period of a supercycle increases from 2''^"^ to 
2^ by a shift the control parameter value from ~ppf_i to 
jLff, the preimage structure advances one stage of compli- 
cation in its hierarchy. Along with this, and in relation to 
the time evolution of the ensemble of trajectories, an ad- 
ditional set of 2^ smaller phase-space gaps develops and 
also a further oscillation takes place in the corresponding 
rate W{t) for finite-period attractors. At /i /ioo the 
time evolution tracks the period-doubling cascade pro- 
gression, and every time t increases from 2^"^ to 2^ the 
fiow of trajectories undergoes equivalent passages across 
stages in the itinerary through the preimage ladder struc- 
ture, in the development of phase-space gaps, and in log- 
arithmic oscillations in W{t). In Fig. 22 we show the 



correspondence between these features quantitatively. 
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FIG. 20; Left panel: Rate W{t) of approach to the Feigen- 
baum attractor a.t fi — fj,oo on logarithmic scales. The curves 
correspond to the values given for the number r of trajecto- 
ries per box, and the expression shown corresponds to power 
law decay with logarithmic oscillations. Right panel: Super- 
position of W{t), for r = 100, on itself via the rescaling shown 
(the same as in Fig. 191 for the horizontal x and vertical y 




FIG. 22: Correspondence between the power law decay with 
log-periodic oscillation features of the rate W{t) with the se- 
quential opening of phase-space gaps. Top panel: The solid 
line is W(t) from Fig. 

" " "■ iV = 



21 and the open circles values are ob 
tained for W{t) from Eq. S at times t = 2^ , N = 1,2, . . 



See text. 
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FIG. 21: Same as in the right panel of Fig. |20| but obtained 
with increased precision. In the left panel the number of ini- 
tial conditions is Nc— 10^, while in the right panel A'^c= 10^. 
The distance A is A = (1 + | - l/a|)/(l -f- | - 1|), where a 
is Feigenbaum's constant. This stems from the fact that all 
initial conditions out of the interval (— l/a, 1) take a value in- 
side this interval in the first iteration. As can be observed, the 
scaling factor for the horizontal axis converges to the exact 
value X — 2. 



IV. g-DEFORMED 
STATISTICAL-MECHANICAL STRUCTURE 

The rate W{t), at the values of time for period dou- 
bling, can be obtained quantitatively from the supercy- 
cle diameters dN,m- Specifically, W{t) = A Zr, A = 



(1 



i)/2 



T 



t — tn, and 



Zr = 



m— 



'N.m 1 



T^2 



N-l 



7V=1,2,3, 



(3) 



Eq. (|3| is an explicit expression equivalent to the nu- 
merical procedure followed in Ref. by the use of 
the triadic cantor set construction of the Feigenbaum at- 
tractor to evaluate the power law exponent tp, and from 
which the value for ip = 0.800138194 is reported. In Fig. 
22 we have added the results of a calculation of W{t) at 
times T = 2^, iV = 1, 2, . . ., according to Eq. (js]). 

The predicted statistical-mechanical structure is ex- 
posed when we identify the (scaled and shifted) decay 
rate Zt as a partition function. From this viewpoint the 
diameters dN,m are configurational terms that in a true 
statistical-mechanical theory would be expected to obey 
some well-defined statistical weights, such as Boltzmann 
factors. To check on this we proceed to determine their 
time dependence given by the bifurcation index N . The 
dN,m scale (asymptotically) with N for fixed as 



N large 



(4) 



where the ay are universal constants obtained, for 
instance, from the finite jump discontinuities of 
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Feigenbaum's trajectory scaling function cr(y) = 
lim„^oo dN.jn+i/dN,m, V = linijv^oo ^2]. The two 

largest discontinuities of a{y) correspond to the thin- 
ner and the fuller regions of the multifractal attrac- 
tor, and for these two regions we have, respectively, 
dNS) ^ and dN,i ^ ^-^(JV-i). [The length of 

the first diameter is di_o = 1 and the equality in Eq. Q 
is approached rapidly with increasing N.] The power 
law in Eq. Q can be rewritten as a g-exponential 
{expij{x) = [1 ~ {q — l)x]~^^^'^~^^) via use of the iden- 
tity A-^+i = {1 + (3 = 2^-1 - 1, that is, 



(5) 



where Qy = I + i^y \ v, 



y — lnQ;j,/ln2, and /3 
1. Similarly, the partition function Z^- ~ t^'^ (or 

2^^^, can be 



2^-1 _ 

Zr ~ e"^+^), with (fi 
expressed as 



T - 1 



In e/ In 2 and r 



Zt ^ expQ(-/?(p), 



(6) 



where Q = 1 + ip^^ and once more /? = r — 1 = 2^^^ — 1. 

Our main contention becomes evident when Eqs. ([5| 
and (|6| are used in Eq. ([s]), to yield 



expg(-/3v3) ~ ^expg^(-/3i/j,). 



(7) 
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FIG. 23: Sector of the bifurcation tree for the logistic map 
(x) that shows the formation of a Pascal triangle of diame- 
ter lengths according to the scaling approximation explained 
in the text, where a ~ 2.5091 is the pertinent universal con- 
stant. 



Eq. ([t]) is similar to a basic statistical-mechanical ex- 
pression; the quantities in it take the following parts: /3 
an inverse temperature, f a free energy (or the prod- 
uct s = —Pf a Massieu thermodynamic potential, or 
entropy), and the Vy configurational energies. However, 
the equality involves g-deformed exponentials in place of 
ordinary exponential functions that would be recovered 
when Q = (/m = 1- It is worth noticing that there is a 
multiplicity of g-indices associated to the configurational 
weights in Eq. Q; however, their values form a well- 
defined family [10 determined by the discontinuities of 
Feigenbaum's function a. 

In the spirit of a mean field approximation we sup- 
pose that, for a given value of N, e.g., N = 3, the di- 
ameters d^.m that are of similar lengths have actually 
equal length, and this length is obtained from those of 
the shortest or longest diameters via a simple scale fac- 
tor; e.g., ds^s = d^^2 = OL~^d^Q = ad-j, i. This initiates 
a degree of degeneracy in the lengths that subsequently 
spreads all the way through the bifurcation tree. See Fig. 



23 [24]. As a result of this approximation, the dAr ,„ scale 
with increasing TV according to a binomial combination of 
the scale factors of those diameters that converge to the 
most crowded and most sparse regions of the multifractal 
attractor. To be precise, the 2^~^ diameters at the A^- 
th supercycle have lengths equal to a~*^^~^~'^a 
occur with multiplicities where A; = 0, 1, .. 



'^^ and 



As shown in Fig. 23 the diameters create a Pascal trian- 
gle across the bifurcation cascade. This feature greatly 
simplifies the evaluation of the partition function and di- 
rectly yields 



Zr 



E 

/c=0 



N - 1 



-2fe 



(a- 



= 2 



We obtain 



(8) 

^ 1.7883, 



e = [ot'^ + a' ^ 
if — 0.8386, and Q — 2.1924, a surprisingly good ap- 
proximation when compared to the numerical estimates 
if = 0.8001and Q = 2.2498 of the exact values. Under 
this approximation all the indices qy in Eq. fzl) are equal, 
qy ^ q ^ \ + , V = lna/ln2, and Eq. (ITrbecomes 



expQ(-/3^) = ^ r!(A^ - 1, k) exp,(-/3i.), 

fc=0 



(9) 



where ^{N — l,k) = a~' {^k^)- thermodynamic lan- 
guage, the approach to the attractor described by Eq. 
([t]) or ([9]) is a cooling process /3 00 in which the free 
energy (or energy) ip is fixed and therefore the entropy 
s = —f3ip is linear in (3. It is instructive to define an 
"energy landscape" for the Feigenbaum attractor as be- 
ing composed of an infinite number of "valleys" whose 
equal-valued minima at (3 00 coincide with the points 
of the attractor on the interval [— a~^,l] [23]. When 
P — 2^^^ — 1, A^ finite, the valleys merge into 2^~^ 
intervals of widths equal to the diameters dN,m- 

As established some time ago, the so-called thermo- 
dynamic formalism is built around the statistical- 
mechanical framework followed by the geometric proper- 
ties of multifractals. The partition function formulated 
to study their properties, like the spectrum of singulari- 
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ties f{a) m, is 

M 

Z(f,q)^^p^/-^ (10) 

m 

where the /,„ (in one-dimensional systems) are M dis- 
joint interval lengths that cover the multifractal set and 
the Pm are probabilities given to these intervals. The 
standard practice consists of demanding that Z(t, q) nei- 
ther vanishes nor diverges in the limit /,„ — for all m 
(and consequently M — > oo). Under this condition the 
exponents r and q define a function -r(q) from which f{a) 
is obtained via Legendre transformation [1 . When the 
multifractal is an attractor its elements are ordered dy- 
namically, and for the Feigenbaum attractor the trajec- 
tory with initial condition xq = generates in succession 
the positions that form the diameters, generating all di- 
ameters djv,m for N fixed between times r — 2^"^ and 
T — 3 X 2^~^. Because the diameters cover the attrac- 
tor it is natural to choose the covering lengths at stage 
N to be /rn''' — dN,in and to assign to each of them the 

same probability ^ = 1/2. For example, within the 
two-scale approximation to the Feigenbaum multifractal 
[T], if^ = a-(^-i-'=)a-2'=, the condition Z(f,q) = 1 re- 
produces Eq. ([sl when pm^ = t^^ = 2^^+^, with t = 1 
and q = — (y3. It should be kept in mind that the "static" 
partition function Z(t, q) is not meant to distinguish be- 
tween chaotic and critical (vanishing A) multifractal at- 
tractors as we do here. As we emphasize below, it is 
the functional form of the link between the probabilities 

pi!^^ and actual time r that determines the nature of the 
statistical mechanical structure of the dynamical system. 

A crossover to q = 1 ordinary statistics when a critical 
attractor turns chaotic; this can be explained as follows. 
We first recall that multifractal sets and their statistical- 
mechanical properties can be retrieved by means of the 
recursive method of backward iteration of chaotic maps 
|25j . A chaotic unimodal map has a two- valued inverse 
and given a position x — Xn binary tree is formed under 
backward iteration, so there are 2" initial conditions xq 
for trajectories that lead to Xn- Since in this case the 
Lyapunov exponent is positive A > 0, lengths expand 
under forward iteration according to / ~ exp(An) and 
contract under backward iteration as / ~ exp(— An). We 
can define, as above, a set of covering lengths Dn^m = 
Amexp(— An), where m relates to the initial condition 
Xq and use of them in a partition function like that in 
Eq. ^ gives 

exp(-/3v5) = ^ A,^exp(-/3A), (11) 

m 

where now (3 — n. Keeping in mind Pesin's theorem Lp is 
plainly identified as the Kolmogorov-Sinai entropy. Now, 
the crossover from g-deformed statistics to ordinary q = 1 
statistics can be observed for control parameter values in 
the vicinity of the Feigenbaum attractor, > /Zoo, when 



the attractor consists of 2" bands, n large. The Lya- 
punov coefficient A of the chaotic attractor decreases with 
= /i— /ioo as A oc 2^" ~ A/i", k — ln2/ InSp, where 
(5i? is the Feigenbaum constant that measures the rate of 
development of the bifurcation tree in control parameter 
space [15]. The chaotic orbit consists of an interband 
periodic motion of period 2" and an intraband chaotic 
motion. The expansion rate '^I^q '^'n\dff^{xi) /dxi\ fluc- 
tuates with increasing amplitude as In r for r < 2" but 
converges to a fixed number that grows linearly with t 
for T ^ 2" [2]. This translates as dynamics with q ^ I 
for r < 2" but ordinary dynamics with g = 1 for r ^ 2". 



V. SUMMARY 

As stated, the dynamics of critical attractors in low- 
dimensional nonlinear maps is a suitable phenomenon for 
assessing the limits of validity and generalizations of or- 
dinary statistical mechanics. There are now positive in- 
dications that the multifractal critical attractors present 
in these maps play this role, since, as it turns out, the 
two sets of dynamical properties — inside and towards the 
Feigenbaum attractor — appear combined in a q-deformed 
statistical-mechanical structure [12]. To obtain this re- 
markable property it is necessary to have access to de- 
tailed information for these two different types of proper- 
ties. The dynamics at the attractor (for both trajectories 
and sensitivity to initial conditions) has been analyzed 
in detail before for period doubling [H [10] and for the 
quasiperiodic [26 transition to chaos. But the dynamics 
on the way to the attractor is only now offered as far as 
we know. 

To begin with, we studied the properties of the first few 
members of the family of superstable attractors of uni- 
modal maps with quadratic maxima and obtained a pre- 
cise understanding of the complex labyrinthine dynamics 
that develops as their period 2^ increases. The study is 
based on the determination of the function tf^xo), the 
time of flight for a trajectory with initial condition xq 
to reach the attractor or repellor. The function tf^xo) 
was determined for all initial conditions xq in a partition 
of the total phase-space —l<xo<l, and this pro- 
vides a complete picture for each attractor-repellor pair. 
We observed how the fractal features of the boundaries 
between the basins of attraction of the positions of the 
periodic orbits develop a structure with hierarchy, and 
how this in turn reflects on the properties of the trajec- 
tories. The set of trajectories produces an ordered flow 
towards the attractor or towards the repellor that reflects 
the ladder structure of the sub-basins that constitute the 
mentioned boundaries. As 2^ increases there is sensitiv- 
ity to the final position for almost all xg, and there is 
a transient exponentially-increasing sensitivity to initial 
conditions for almost all xq. We observed that transient 
chaos is the manifestation of the trajectories' controlled 
flow out of the fractal boundaries, which suggests that for 
large 2^ the flow becomes an approximately self-similar 
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sequence of stages. As a final point, in the Appendix, 
we look at the closing segment of trajectories at which a 
very fast convergence to the attractor positions occurs. 
We found "universality class" features, as the trajecto- 
ries and sensitivity to initial conditions are replicated by 
a RG fixed-point map obtained under functional compo- 
sition and rescaling. This map has the same q-deformed 
exponential closed form found to hold also for the pitch- 
fork and tangent bifurcations of unimodal maps (THl HZ] . 

Subsequently, we examined the process followed by an 
ensemble of uniformly distributed initial conditions Xq 
across the phase-space to arrive at the Feigenbaum at- 
tractor, or get captured by its corresponding repellor. 
Significantly, we gained understanding concerning the 
dynamical ordering in xq, in relation to the construc- 
tion of the families of phase-space gaps that support the 
attractor and repellor, and about the rate of approach of 
trajectories towards these multifractal sets, as measured 
by the fraction of bins W{t) still occupied by trajecto- 
ries at time t. An important factor in obtaining this 
knowledge has been the consideration of the equivalent 
dynamical properties for the supercycles of small periods 
in the bifurcation cascade. As we have seen, a doubling 
of the period introduces well-defined additional elements 
in the hierarchy of the preimage structure, in the family 
of phase-space gaps, and in the log-periodic power law 
decay of the rate W{t). We have then corroborated the 
wide-ranging correlation between time evolution at /ioo 
from t — wp to t —^ oo with the static period-doubling 
cascade progression from /i = up to /loo • As a result of 
this we have acquired an objective insight into the com- 
plex dynamical phenomena that fix the decay rate W{t). 
We have clarified the genuine mechanism by means of 
which the discrete scale invariance implied by the log- 



periodic property in W{t) arises, that is, we have seen 
how its self-similarity originates in the infinite hierarchy 
formed by the preimage structure of the attractor and re- 
pellor. The rate W{t) can be obtained quantitatively [see 
Eq. ([3|] from the supercycle diameters dM,m- These ba- 
sic data descriptive of the period-doubling route to chaos 
are also a sufficient ingredient in the determination of 
the anomalous sensitivity to initial conditions for the dy- 
namics inside the Feigenbaum attractor |10j . 

Finally, the case is made that there is a statistical- 
mechanical property underlying the dynamics of an 
ensemble of trajectories en route to the Feigenbaum at- 
tractor (and repellor). Eq. ^ is identified as a partition 
function built of g-exponential weighted configurations, 
and in turn, the fraction Zr of phase-space still occupied 
at time r is seen to have the form of the g-exponential 
of a thermodynamic potential function. This is argued 
to be a concrete, clear and genuine manifestation of 
q-deformation of ordinary statistical mechanics where 
arguments can be made explicit and rigorous. There is 
a close resemblance with the thermodynamic formalism 
for multifractal sets, but it should be stressed that 
the deviation from the usual exponential statistics is 
dynamical in origin, and due to the vanishing of the 
(only) Lyapunov exponent. 
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VI. APPENDIX: SUPER-STRONG 
INSENSITIVITY TO INITIAL CONDITIONS 

The ordinary Lyapunov exponent A 



limt^oo t ^ In dfj^^ {x)/dx 



for the supercycle at- 

tractors diverges to minus infinity [dfj^^{x)/dx = at 
X — 0], therefore the sensitivity to initial conditions 
cannot have an exponential form = exp(At) with 
A < 0. This appendix is a brief account of the deter- 
mination of at the closing point of approach to the 
supercycle attractors. 




a numerical investigation are shown in Fig. 24 Only the 
first two steps of a trajec tor y with xq = 0.1 of the map 



(2 ) 

' can be seen in Fig. 



24 a). A considerable enlarge- 
ment of the spatial scale (which requires computations 
of extreme precision |27|) makes it possible to observe a 
total of seven steps (56 iterations in the original map), as 
shown with the help of logarithmic scales in Fig. |24j^b). 
Fig. 24 ^c) shows both the same trajectory and the sen- 
sitivity = dxt/dxo in a logarithmic scale for Xt and 
ft and a normal scale for the time t. The trajectory is 
accurat ely reproduced [indistinguishable from the curve 
(c)] by the expression xt — exp{bexpct), 

In 2, where u > is 



24 



in Fig 

b = Inuxo (with Xq > 0), and c 



obtained from the form | /- 



(2") I 



taken by the 2^-th 



composed map close to a; = 0. This expression for Xt is 
just another form of writing uxt = (uxq)'^ , the result 
of repeated iteration of ux'^. For the logarithm of the 
sensitivity we have In ft = — In + i In 2 + In Xt where 
the last (large negative) term dominates the first two. 
Thus, we find that the sensitivity decreases more quickly 
than an exponential, and, more precisely, decreases as 
the exponential of an exponential. 

We note that, associated with the general form 

fjj (x) ~ ux of the map in the neighborhood of a; = 0, 
there is a map f*{x) that satisfies the functional compo- 
sition and rescaling equation f*{f*{x)) — a~^f*{ax) for 
some finite value of a and such that /* (x) ~ ux^ + o{x'^). 
The fixed-point map /* (x) possesses properties common 
to all superstable attractors of unimodal maps with a 
quadratic extremum. Indeed, there is a closed form ex- 
pression that satisfies these conditions, which is f*{x) = 
a;expg(u'~^2;). The fixed-point map equation is satisfied 
with a = 2i/(9"i) and q= 1/2. The same type of RG 
solution has been previously found to exist for the tan- 
gent and pitchfork bifurcations of unimodal maps with 
general nonlinearity z > 1 [TH [171 IE] ■ Use of the map 
f*{x) reproduces the trajectory Xi/xq and sensitivity ft 
shown in Fig. 24 (c). 



FIG. 24: Panel (a): Detail of a trajectory in its final stage of 
approach to the x = Q attractor position (dotted line) for the 

map /p {x) (solid line). In this example N — 3. Panel (b): 
Same as (a) in double logarithmic scale. Panel (c): Trajectory 
Xt and sensitivity to initial conditions ft in logarithmic scale 
versus time t. Both functions are indistinguishable. Note that 
there is an ultra-rapid convergence, of only a few time steps, 
to the origin a; = 0. See text. 



Representative results for the last segment of a trajec- 
tory and the corresponding sensitivity ft obtained from 



